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To examine whether microbial community structure differs across rice genotypes, automated ribosomal intergenic 
spacer analysis (ARISA) was conducted. Nine cultivars of Oryza saliva ssp. indica or japonica and seven lines of 
other Oryza species were grown in paddy fields with low, standard, and high levels of N fertilization. Multidimensional 
scaling plots of bacterial ARISA for aerial parts of rice (shoots) revealed that the structure of shoot bacterial communities 
was significantly affected by plant genotype {indica or japonica) based on similarity tests, whereas root microbiomes 
were largely affected by the N fertilization level. 
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MICROBES AND 
ENVIRONMENTS 



A wide range of microorganisms, including bacteria and 
fungi, have been found in the phyllosphere and rhizosphere 
(21, 25, 26), and reside in and on plants as endophytes 
and epiphytes (22, 25, 30). These symbiotic microbes could 
be important components of the proximate mechanisms 
underlying plant functional traits such as nutrient acquisition, 
plant defense, plant morphology and abiotic stress tolerance 
(7); however, many questions remain about the driving forces 
for shaping a community structure of plant-associated 
microbes (9). 

In studies of how plant genetic factors control microbial 
community structure, artificially generated genetic variations 
are an important basic resource (11). In rice, Ikeda et al. (14) 
reported significant impacts of the OsCCaMK gene on the 
diversity of root-associated bacteria under both paddy and 
upland field conditions using OsCCaMK mutants screened 
from a Tosl7 mutant panel (2). CCaMK plays an important 
role in a common symbiosis pathway that leads to successful 
rhizobial and arbuscular mycorrhizal symbioses in plants 
(18). In addition, naturally occurring genetic variation in 
plants is another important resource (17). Genotypes derived 
from natural plant variations have been shown to affect root- 
associated bacterial communities in rice (10). Hardoim et al. 
(10) demonstrated that plant genotypes helped to shape 
bacterial community structures in rice roots; however, the 
effects of plant genotype on aerial parts of rice (shoot)- 
associated bacterial communities have not been analyzed by 
culture-independent methods. 

Nitrogen is the most important mineral nutrient for crop 
production, and an adequate supply of nitrogen fertilizer is 



* Corresponding author. E-mail: kiwamu@ige.tohoku.ac.jp; 

Tel: +81-22-217-5684; Fax: +81-22-217-5684. 
f Present address: Plant Breeding, Genetics, and Biotechnology, 

International Rice Research Institute, DAPO Box 7777, Metro 

Manila, The Philippines 



essential for sustaining high yields. Previous studies have 
shown that nitrogen application influences rice-associated 
microbial communities using polymorphisms of 16S rRNA 
and nifH genes (5, 15, 29). 

In general, genetic fingerprinting techniques allow high 
throughput and comparative profiling of numerous samples 
for microbial community analyses (3, 26). Automated 
ribosomal intergenic spacer analysis (ARISA) is a highly 
sensitive method due to the laser detection of fiuorescently 
labeled DNA, and it is a rapid and effective technique that 
can be used in conjunction with more accurate but labor- 
intensive methods (26). Using this method, Ikeda et al. (13) 
showed that the bacterial community structure in soybean 
roots could be classified into three groups according to the 
plant genotype. 

In the present study, ARISA was used to reveal the effects 
of inter- and intraspecific genetic variations in the genus 
Oryza on bacterial communities. Nine cultivars of O. sativa 
and seven lines of other Oryza species were used for shoot- 
associated bacterial community analysis, and five cultivars 
of O. sativa were used for root-associated bacterial commu- 
nity analysis (Table I). 

A total of 48 seedlings of each cultivar and line were 
planted and grown in paddy fields at Kashimadai, Miyagi, 
Japan, as described by Obara et al. (24). Low nitrogen (LN) 
and standard nitrogen (SN) fields were managed with 0 and 
30 kg N ha ' of nitrogen fertilization, respectively. In the 
high nitrogen (HN) field, 270 kg N ha ' of ammonium sulfate 
(Ube Agri-Materials, Tokyo, Japan) was applied in addition to 
basal fertilizer. To reach 270 kg N ha ', 30 or 60 kg N ha ' 
of ammonium sulfate was applied as additional fertilization 
every 2 weeks. Plants were sampled 90 days after trans- 
planting. Each rice plant was dug up in a square (30 cm x 30 
cm) to a depth of about 30 cm from the soil surface. The 
roots were washed with running tap water in the laboratory 
until the soil was removed from roots, and around 20 g shoots 
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Table 1. Rice plant materials used in this study 



Cultivar/line 


Species 


Genome^ 


Nipponbare 


O. SQtiva L. ssp. japonicu 


AA 


Sasanishiki 


O. sativa L. ssp. japonico 


AA 


Taichung 65 


O. sativa L. ssp. japonica 


AA 


Gemdjah Beton 


O. sativa L. ss^. japonica 


AA 


Koshihikari 


O. sativa L. ^s^. japonica 


AA 


Habataki 


O. sativa L. ssp. indica 


AA 


Kasalath 


O. sativa L. ssp. indica 


AA 


1R24 


O. sativa L. ssp. indica 


AA 


1R36 


O. sativa L. ssp. indica 


AA 


WK21 


O. glaberrima Steud. 


AA 


W106 


O. ruflpogon sensu lato 


AA 


W1965 


O. ruflpogon sensu lato 


AA 


W1967 


O. ruflpogon sensu lato 


AA 


W630 


O. ruflpogon sensu lato 


AA 


W1515 


O. punctata Kotschy ex Steud. 


BB 


W1527 


O. eichingeri Peter 


CC 



" Cultivars and lines were classified based on the genome composition 
according to the degree of meiotic chromosome pairing in hybrid plants 
(23). 



and a whole root in each plant were stored separately at 
-80°C until molecular analysis. A total of 20 g shoots and 
whole root tissue from one plant were separately ground to 
a powder in liquid nitrogen with a mortar and pestle. DNA 
was extracted from 200 to 300 mg powdered tissue by the 



DNA extraction method developed by Ikeda et al. (12). The 
remaining shoot samples were dried and then digested with 
sulfiiric acid to quantify total nitrogen by the Kjeldahl method 
(Table SI). 

The multidimensional scaling (MDS) method was used in 
conjunction with ARISA to evaluate similarities of bacterial 
communities in shoot and roots. The bacterial primer set 
ITSF/ITSReub (3) was used for ARISA. PCR amplification 
was carried out as described by Ikeda et al. (13). PCR products 
(1 ^iL) were mixed with 1 \\L LIZ 1200 internal size standard 
(Applied Biosystems, Foster City, CA, USA), and 20 |aL 
deionized formamide was added. The mixture was denatured 
at 95°C for 5 min and cooled on ice. Next, the PCR product 
was placed in an ABI 3730x1 DNA Analyzer (Applied 
Biosystems). The profile data (up to 1,200 bases) obtained 
by ARISA were initially analyzed with ABI GeneMapper 
software (Applied Biosystems) and were processed further 
with the RIBOSORT program (27) to assign fragment size 
and calculate the relative abundance of each ribotype. 
Fragments with fewer than 100 fluorescence units were 
eliminated from the analyses. Using the R program (27), an 
MDS plot was generated from the ribotypes using a similarity 
matrix produced by the RIBOSORT program with default 
parameters and the vegan package with the Bray-Curtis index 
of dissimilarify. Analyses of similarities (ANOSIM) were 
performed to test for significant differences using 1000 
permutation tests. The resulting test statistic R indicates the 
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Fig. 1. Multidimensional scaling plots generated from ARISA profiles with primers ITSF and ITSReub for rice (A) shoot- and (B) root- 
associated bacterial communities under low (LN), standard (SN), and high (HN) nitrogen conditions. Nine cultivars of O. sativa and seven lines of 
other Oryza species were used for shoot-associated bacterial community analysis with three replications, and five cultivars of O. sativa were used 
for root-associated bacterial community analysis with at least two repUcations (Table 1). 
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Fig. 2. Multidimensional scaling plots generated from ARISA profiles with primers ITSF and ITSReub for rice shoot-associated bacterial com- 
munities under (A) low, (B) standard, and (C) high nitrogen conditions. _/(3po«/ca and indica indicate O. sativa spp. japonica and indica, respec- 
tively (Table 1). W106, W1967, W1965 and W630 belongs to O. mflpogon (Table 1) («=3). 



degree of separation, ranging from 0 (no separation) to 1 
(complete separation). Numbers of operational taxonomic 
units in each sample are shown in Table S2. 

MDS was performed to compare the complex ARISA 
profiles of the bacterial community structures among the 
plant genotypes within the genus Oryza (Table 1). In O. 
sativa, MDS plots of the bacterial ARISA of shoot samples 
indicate that the bacterial communities were apparently 
separated between japonica and indica groups (Fig. lA), 
which were statistically significant by similarity tests (ANO- 
SIM) (Table S3). MDS plots of the bacterial ARISA of 
shoot samples from other Oryza species {O. glaberrima, O. 



rufipogon, O. punctata, and O. eichingeri) (Table 1) did not 
form a distinct cluster, but were widely dispersed (Fig. lA), 
suggesting that the diversity of bacterial community in shoots 
of O. sativa was lower than those of other Oryza species 
(Fig. lA). This may reflect the presence of more genetic 
variations found in these Oryza species than in those of O. 
sativa (16). 

MDS plots of the bacterial ARISA of shoot samples 
indicated that the community structures for japonica and 
indica genotypes formed clear clusters along dimension 1 

under all nitrogen conditions (Fig. 2), suggesting that the 
effect of a plant genotype in O. sativa species on the shoot- 
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associated bacterial community is greater than that of the 
nitrogen fertilization level. Although the nitrogen concentra- 
tions of rice shoots under HN conditions were apparently 
higher than those under other conditions (SN and LN), no 
correlation was observed between the N concentration and 
microbial community (Table SI). 

Diverse environmental factors control the establishment of 
microbial communities in the phyllosphere, but recently it 
has been recognized that a plant genotype plays an important 
role in selecting phyllosphere communities (30). There are 
considerable variations between indica and japonica cultivars 
in the concentrations of chemical components, such as nitrate 
and heavy metals (1,6). In addition, De Costa et al. (4) found 
that leaves from different rice genotypes have different 
culturable bacterial communities and concluded that the 
difference was significantly correlated with anatomical and 
physiological differences, such as leaf hair length and density, 
leaf temperature, stomatal density, and transpiration rate. 
Among these ttaits, the variation between indica and japonica 
cultivars in stomatal density and size and canopy temperature 
is well recognized, and the genetic factors underlying these 
traits have been identified (19, 28). Thus, these variations in 
the chemical characteristics and microstructures between 
indica and japonica might affect shoot-associated bacterial 
communities. 

For root-associated bacterial community analysis, Sasa- 
nishiki, Taichung 65, Gemdjah Beton, IR24, and IR36 were 
examined (Table S2). MDS plots for bacterial ARISA of root 
samples indicated that root bacterial communities under SN 
and HN conditions formed a tight cluster (Fig. IB), suggesting 
that the effect of nitrogen fertilization on the bacterial 
community associated with roots was greater than that of the 
plant genotypes examined. Demba Diallo et al. (5) showed 
that nitrogen treatment had a strong effect on the composition 
and diversity of expressed nifH pools that shifted towards 
methylotroph-related nitrogenases. In another study, coloni- 
zation of sugarcane by Acetobacter diazotrophicus has been 
inhibited by high nitrogen fertilization (8). The bacterial 
community of bulk soil would help to understand the effect 
of nitrogen fertilization on the paddy field ecosystem (31). 

Different genotypes of Arabidopsis thaliana have been 
used for structural comparisons of root microbial communi- 
ties to improve our imderstanding of plant-microbe interac- 
tions (20). To the best of our knowledge, the present study 
is the first report that shoot-associated microbial communities 
are largely dependent on rice genotypes (O. sativa ssp. 
japonica and indica). 
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